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Abstract Sulfur is an abundant element which remains 
undepleted in diffuse interstellar gas (Ay < 1) but it 
is traditionally assumed to deplete on dust grains at 
higher densities and larger Ay. 

Photodissociation regions (PDRs) are an interest- 
ing intermediate medium between translucent and dark 
clouds where the chemistry and energetics are domi- 
nated by the illuminating FUV radiation field. Thus 
they can provide new insights about the sulfur deple- 
tion problem. However, physical and chemical gradients 
in PDRs take place at small angular scales (~1 to 10"). 
Aperture synthesis observations are therefore required 
to resolve such gradients. Besides, a complete under- 
standing of molecular excitation is needed to correctly 
determine molecular abundances but also the preavail- 
ing physical conditions. Hence, multi-J observations at 
increasing frequencies are also required. Such high an- 
gular resolution and broad frequency coverage observa- 
tions will be provided by ALMA in the near future. 

In this work we present IRAM-PdBI observations of 
the CS J=2— 1 line toward the Horsehead PDR comple- 
mented with IRAM-30m observations of several rota- 
tional lines of different sulfur bearing molecules (CS, 
HCS + , SO, H2S, etc.). Photochemical and nonlocal, 
non-LTE radiative transfer models adapted to the Horse- 
head geometry have been developed. The gas phase sul- 
fur abundance has been inferred in the PDR. 
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1 Introduction 

Sulfur is an abundant element (the solar photosphere 
abundance is S/H=1.38x 10 -5 ), which remains unde- 
pleted in diffuse interstellar gas and H11 regions but it 
is traditionally assumed to deplete on grains in higher 
density molecular clouds by factors as large as ^1000. 
This conclusion is simply reached by adding up the ob- 
served gas phase abundances of S-bearing molecules in 
well known dark clouds such as TMC1. As sulfur is eas- 
ily ionized (ionization potential ^10.36 eV), sulfur ions 
are probably the dominant gas-phase sulfur species in 
translucent gas. Ruffle et al. (1999) proposed that if 
dust grains are typically negatively charged, S + may 
freeze-out onto dust grains during cloud collapse more 
efficiently than neutral species such as oxygen. How- 
ever, the nature of sulphur on dust grains (either in 
mantles or cores) is not obvious. The absence of strong 
IR features due to S-bearing ices in many ISO's mid- 
IR spectra and the presence of S 11 recombination lines 
in dark clouds such as Rho Ophiuchi all argue against 
a large depletion of sulfur from the gas phase. Besides, 
the abundance of species such as CS or H 2 S may indi- 
cate that something important is lacking from chemi- 
cal models or that an abundant sulfur-bearing carrier 
has been missed. Therefore, the abundances of sulfur 
species remain interesting puzzles for interstellar chem- 
istry. PDRs offer an ideal intermediate medium between 
diffuse and dark cloud gas to investigate the sulfur de- 
pletion problem. In particular, the Horsehead western 
edge is a PDR viewed nearly edge-on. The intensity of 
the incident FUV radiation field is x — 60 relative to 
the interstellar radiation field (ISRF) in Draine's units. 
At a distance of 400 pc, the Horsehead PDR is very 
close to the prototypical source needed to serve as a 
reference to models (Pety et al. 2006). 
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Since 2001 we have been studying the Horschcad 
nebula at high angular resolution (< 5") with the IRAM 
PdBI interferometer, and complemented with larger scale 
and higher frequency observations with the IRAM-30m 
single-dish telescope. Following these and other works, 
the Horsehead PDR is now considered a well under- 
stood source, ideal to analyse the chemical stratification 
predicted in clouds illuminated by a FUV radiation field 
(see Pety et al. these proceedings). These studies have 
led e.g., to a significant increase of our knowledge of the 
carbon chemistry (from PAHs to small hydrocarbons) 
in FUV irradiated gas (see Gerin et al. 2005; Pety et 
al. 2005). More recently, we have tried to accurately 
determine and understand the observed CS and HCS + 
abundances in the Horschcad PDR as a tool for estimat- 
ing the sulfur gas phase abundance. Compared to other 
works, we now specifically treat the nonlocal, non-LTE 
molecular excitation and radiative transfer along the 
cloud edge using the output predictions of photochem- 
ical models. In this way, we can consistently link the 
thermal and chemical predictions of PDR models with 
the absolute molecular line intensities detected by tele- 
scopes. This, sometimes difficult, step has a crucial im- 
portance in the accurate determination of abundances 
and in our understanding of the prevailing chemistry 
and physical conditions by a direct comparison of PDR 
models with molecular line emissivities. 



In this work we present 3.65" x 3.34" angular reso- 
lution IRAM Plateau de Bure Interferometer observa- 
tions of the CS J— 2-1 line (Fig. 2 left), complemented 
with IRAM-30m observations of several rotational lines 
of CS, C 34 S, HCS+, SO, 34 SO, H 2 S, etc. (some of them 
shown in Fig. 1). Beam sizes of single-dish observations 
range from ~10" at 1 mm to ~30" at 3 mm. Data re- 
duction was done with the IRAM GILD AS software. 



2 Numerical methodology 

To analyse the gas phase sulfur chemistry we use the 
Meudon PDR code, a photochemical model of a unidi- 
mensional stationary PDR (Le Petit et al. 2006). In few 
words, the PDR code solves the FUV radiative trans- 
fer in an absorbing and diffusing medium of gas and 
dust (Goicoechea & Le Bourlot 2007). This allows the 
explicit computation of the FUV radiation field (con- 
tinuum+lines) and therefore, the explicit integration 
of consistent C and S photoionization rates together 
with H 2 , CO, 13 CO, and C 18 photodissociation rates. 
Our standard conditions for the model of the Horsehead 
PDR include a power-law density profile and a FUV ra- 
diation field enhanced by a factor x = 60 with respect 
to the Draine ISRF. Different sulfur gas phase abun- 
dances, S/H, have been investigated (see Fig. 2 right). 
To be consistent with PdBI CO multi-line observations, 
the thermal balance was solved until the gas tempera- 
ture reached a minimum value of 30 K, then a constant 
temperature was assumed. 

For the molecular excitation and radiative transfer 
in the PDR we use a nonlocal, non-LTE Monte Carlo 
code that allows us to directly compare our millime- 
ter line observations. The code handles spherical and 
plane-parallel geometries and accounts for line trap- 
ping, collisional excitation, and radiative excitation by 
absorption of microwave cosmic background and dust 
continuum photons. Arbitrary density, temperature or 
abundance profiles, and complex velocity gradients can 
be included. A detailed description of the model is given 
in Goicoechea (2003) and Goicoechea et al. (2006). 

The following methodology was carried out: a full 
PDR model with Horsehead standard conditions was 
run with a particular choice of the density gradient. 
Afterwards, the PDR output was used as input for 
the excitation and radiative transfer calculation. In this 
way, physical parameters can be tuned more accurately 
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Fig. 2 Left panel: CS J=2-l integrated emission map obtained with the PdBI. The map center has been set to the mosaic phase 
center: RA(2000) = 05h40m54.27s, Dec(2000) = -02°28'00". The map size is 110" X 110". The synthesized beam is plotted in the 
bottom left corner. The map has been rotated by 14° counter-clockwise around the image center to bring the exciting star direction 
in the horizontal direction as this eases the comparison of the PDR models. Right panel: Photochemical model predictions for the 
physical and FUV illuminating conditions prevailing in the Horsehead PDR showing the CS and HCS + abundances as a function of the 
sulfur gas phase abundance. Horizontal shaded regions show the CS and HCS + abundances derived from the single-dish observations 
and radiative transfer modeling. Note that for clarity HCS + abundances have been multiplied by a factor of 1000. The shaded vertical 
region shows the estimated sulfur abundance in the Horsehead nebula derived from the constrained fits of CS and HCS + abundances. 



by iteration of different radiative transfer models. Once 
better parameters have been found, a new PDR com- 
putation is performed with this choice of physical pa- 
rameters. We have so far analyzed CS, C 34 S, C 18 and 
HCS + species. Synthetic abundance profiles are consis- 
tently computed as a function of the edge distance Sx 
(in arcsec) and directly compared with interferometric 
observations (see Fig. 3). We are now analyzing a larger 
collection of sulfur bearing molecules. 

3 Main Results 

We have analyzed interferometric CS J =2-1 line maps 
of the Horsehead PDR at a 3.65" x 3.34" resolution 
together with single-dish observations of several rota- 
tional lines of CS, C 34 S and HCS+. We have studied 
the CS photochemistry, excitation and radiative trans- 
fer using the latest HCS + and OCS + dissociative re- 
combination rates (Montaigne et al. 2005) and CS colli- 
sional cross-sections (Lique et al. 2006). The main con- 
clusions are as follows (the whole work can be found in 
Goicoechea et al. (2006)): 

1. CS and C 34 S rotational line emission reveals mean 
densities around 7i(i/ 2 )=(0. 5-1.0) x 10 5 cm~ 3 . The 
CS J=5-4 lines show narrower line widths than 
the low-J CS lines and require higher density gas 
components, ~(2-6)xl0 5 cm~ 3 , not resolved by a 
~10" beam. These values are larger than previous 



estimates based on CO single-dish observations. It 
is likely that dumpiness at scales below ~0.01 pc 
and/or a low density envelope play a role in the CS 
line profile formation. 

2. Nonlocal, non-LTE radiative transfer models of op- 
tically thick CS lines and optically thin C 34 S lines 
provide an accurate determination of the CS abun- 
dance, x(CS)=(7±3)xl0~ 9 . We show that radiative 
transfer and opacity effects play a role in the result- 
ing CS line profiles but not in C 34 S lines. Assuming 
the same physical conditions for the HCS + molecu- 
lar ion, we find x(HCS+) = (4±2) x 10" 11 . 

3. According to photochemical models, the gas phase 
sulfur abundance required to reproduce these CS 
andHCS+ abundances is S/H=(3.5±1.5) x 10~ 6 , only 
a factor ~4 less abundant than the solar elemental 
abundance. Larger sulfur abundances are possible 
if the gas is significantly warmer. Thus, the sulfur 
abundance in the PDR is very close to the undc- 
pleted value observed in the diffuse ISM. The pre- 
dicted CS/HCS+ abundance ratio is close to the ob- 
served value of ~175, especially if predicted HCS + 
peak abundances are considered. If not, the HCS + 
production is underestimated unless the gas is in a 
higher ionization phase, e.g. if the cosmic ray ioniza- 
tion rate is increased by ~5. We are currently work- 
ing on a larger inventory of sulfur-bearing molecules 
detected by us in the Horsehead PDR to test other 
limitations of the sulfur chemical network. 
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Fig. 3 PdBI C ls O J=2— 1 and CS J=2-l spectra along the exciting star direction, Sy = 0" (upper panel) and 8y = 30" (lower panel). 
Monte Carlo radiative transfer models using the output of PDR models for C ls O (blue curve) and CS (red curve) for the physical 
conditions discussed in the text (assuming that the PDR is inclined relative to the line of sight by a ip=5° angle). Modeled line profiles 
have been convolved with an appropriate gaussian beam corresponding to each synthesized beam. 



4. High angular resolution PdBI maps reveal that the 
CS emission does not follow the same morphology 
shown by the small hydrocarbons emission the PDR 
edge (Pety et al. 2005). In combination with previ- 
ous PdBI C 18 observations we have modeled the 
PDR edge and confirmed that a steep density gradi- 
ent (see Habart et al. 2005) is needed to reproduce 
CS and C 18 observations. The resulting density 
profile qualitatively agrees to that predicted in nu- 
merical simulations of a shock front compressing the 
PDR edge to high densities, n.(i?2)— 10 5 cm -3 , and 
high thermal pressures, ~(5-10)xl0 6 K cm~ 3 . 

5. Conventional PDR heating and cooling mechanisms 
fail to reproduce the temperature of the warm gas 
observed in the region by at least a factor ~2. Ad- 
ditional mechanical heating mechanisms associated 
with the gas dynamics may be needed to account 
for the warm gas. The thermal structure of the PDR 
edge is still not fully constrained from observations. 
This fact adds uncertainty to the abundances pre- 
dicted by photochemical models. 

Observational studies of ISM clouds show that many 
physical and chemical variations occur at very small 
angular scales. For PDRs, an accurate molecular in- 
ventory, as a function of the distance from the illumi- 
nating source, can only be obtained from interferomet- 
ric observations. Aperture synthesis observations con- 
tain detailed information about density, temperature, 
abundance and structure of the cloud, but only spe- 
cific radiative transfer and photochemical models for 
each given source are able to extract the information. 
A minimum description of the source geometry is usu- 
ally needed. Future observations with ALMA will allow 
us to finally characterize PDRs at the appropriate spa- 
tial scales of their physical and chemical gradients. 
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